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摘要 : 对 地 平 式 望远镜 机 架 模 型 拟 合 的 指向 误差 数据 展开 分 析 。 介绍 了 指向 误差 数据 的 
获取 方法 。 采 用 最 小 二 乘法 求解 修正 系数 结合 方位 轴 位 置 和 高 度 轴 位 置 拟 合 出 指向 修正 量 。 
在 显著 性 水 平 a=0.05 下 进行 残 差分 析 ， 吻 除 偏离 过 大 的 点 以 提高 拟 合 精度 。 通 过 计算 均 方 
根 误差 ， 相 关系 数 ,F 检验 的 p 值 验证 了 机 架 模 型 有 效 性 。 本 文 说 明了 增加 测试 数据 后 ,未 
经 残 差 检验 不 一 定 能 提高 拟 合 精度 ; 以 及 分 次 获取 的 误差 数据 可 合并 处 理 。 将 求解 的 系数 写 


入 配置 文件 ,在 指向 与 跟踪 过 程 中 根据 光栅 编码 器 的 反馈 值 实时 修正 。 
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Screening and analysis of pointing error data under 
frame model 
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Abstract: The pointing error data of the frame model which apply to the altazimuth telescope 
are analyzed. The method of obtaining pointing error data is introduced in this paper. The least 
square method is used to solve the correction coefficient, and the pointing correction amount is 
fitted by combining the position of the azimuth axis and the height axis. At the significance level 
ofu= 0.05, the residual analysis ls carried out to eliminate those points with excessive deviation 
to improve the fitting accuracy. The validity of the model is verified by calculating RMS error, 
correlation coefficient and p value of F test. This article explains that the fitting accuracy can not 
be improved without residual test after adding test data, and the error data obtained in different 
times can be combined for processing. Coefficients are written into the configuration file, and then 


correct pointing error in real time according to the feedback value of the grating encoder during 
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the pointing and tracking process. 
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了 中 


1 引 


我 国有 许多 学 者 在 天 文 望 远 镜 的 指向 误差 修正 方面 进行 过 研究 ,中科院 南京 天 文 仪器 公 
司 的 朱 庆 生 、 陈 伟 民 等 人 建立 了 机 架 模 型 “， 经 修正 后 的 指向 精度 优 于 5 。 中 国 科学 院 紫 金 
山 天 文 台 的 张 晓 祥 、 吴 连 大 对 静态 指向 模型 的 基本 参数 进行 了 讨论 ， 通 过 逐步 回归 分 析 ， 确 
定 了 码 盘 零 点 差 、 度 盘 偏 心 误差 等 误差 源 作为 显著 回归 因子 ”“。 中 国 科学 院 国家 天 文 台 的 李 
振 伟 、 杨 文 波 等 人 对 比 了 球 谐 函 数 模型 和 水 平 式 望远镜 指向 模型 , 并 通过 水 平 式 望远镜 观测 
验证 了 模型 的 有 效 性 ”。 在 射电 望远镜 的 指向 性 能 改进 上 ， 中 国 科 学 院 新 疆 天 文 台 的 温 浩 兴 
等 人 通过 修正 轨道 高 差 造 成 的 指向 偏差 ， 提高 了 天 线 整体 指向 精度 “。 本 文通 过 软件 的 方式 
修正 指向 误差 , 对 机 架 模型 下 误差 数据 的 获取 及 系数 求解 做 了 说 明 , 并 通过 残 差 分 析 的 方法 
筛选 指向 误差 数据 ,计算 均 方 根 误差 以 验证 模型 对 优化 后 数据 的 内 符合 精度 , 排除 找 星 错误 
等 误 操 作对 修正 结果 的 影响 , 提高 拟 合 精度 。 在 此 基础 上 分 析 加 密 测 量 点 数 对 指向 修正 精度 
的 影响 , 并 说 明了 分 次 获取 的 误差 数据 ， 残 差 检验 后 可 合并 处 理 。 使 用 本 文 的 方法 对 工 米 人 
卫 激 光 测 距 望 远 镜 和 地 基 激 光 发 射 望 远 镜 获 取 的 指向 误差 数据 进行 分 析 , 修正 精度 均 得 到 有 


2 指向 误差 影响 因素 分 析 


天 文 望远镜 在 设计 、 加 工 、 制 造 、 装 配 等 环节 ， 由 于 多 种 因素 引入 指向 误差 ， 导 致 望 远 
镜 主 光 轴 显示 到 达 的 位 置 与 目标 星 实际 的 位 置 存在 偏差 , 在 卫星 跟踪 中 不 能 准确 指向 卫星 的 
预报 位 置 。 这 些 因素 主要 包含 机 架 的 制造 与 安装 误差 、 轴 系 倾斜 、 轴 间 的 不 正 交 、 重 力 变形 、 
镜 简 与 光学 系统 的 弯 沉 “、 光 栅 编 码 器 误差 、 视 准 轴 误 差 等 。 指 向 误差 一 般 用 均 方 根 误差 
(RMSE) 表示 。 

修正 指向 误差 可 采用 硬件 修正 和 软件 修正 两 种 方法 。 便 件 修 正 指 机 械 和 光学 系统 设计 加 
工 、 闭 调 过 程 中 ， 采 用 轴 系 校准 、 视 轴 校 准 、 光 栅 编 码 器 测 角 组 件 校准 等 方法 减 小 误差 。 但 
望远镜 装 调 结束 后 就 难以 通过 硬件 修正 减 小 误差 , 且 随 着 运动 部 件 的 磨损 ， 硬 件 修 正 效果 大 
打折 扣 ， 便 件 修正 效果 有 限 且 费用 较 高 。 软 件 修正 提取 了 影响 望远镜 指向 误差 的 各 个 因素 ， 
利用 系统 误差 重复 性 的 特征 ， 建 立 数 学 模型 并 通过 计算 机 软件 修正 。 
通过 对 望远镜 机 械 系统 的 分 析 , 可 以 估计 出 不 加 修正 前 的 误差 量 , 这 部 分 误差 为 系统 误 
差 。 就 地 平 式 望 远 镜 而 言 ， 指 向 误差 中 系统 误差 较 大 ， 随 机 误差 较 小 ， 由 于 系统 误差 与 望 远 
镜 轴 系 的 角度 有 关 ， 可 建立 高 度 误差 AH 和 方位 误差 44 的 拟 合 模型 。 

望远镜 在 指向 与 跟踪 前 需要 建立 指向 误差 模型 ,通过 回归 分 析 确 定 拟 合 参数 。 在 运行 时 ， 
计算 机 实时 计算 望远镜 在 当前 指向 处 的 系统 误差 估计 值 , 并 对 光栅 编码 器 反馈 的 位 置 值 进 行 
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修改 ， 则 可 消除 指向 误差 中 的 系统 误差 ， 大幅 提高 望远镜 的 指向 精度 和 盲 测 的 跟踪 精度 ， 美 


的 SLR2000 经 模型 修正 后 指向 误差 可 小 于 1”。 


望远镜 指向 误差 修正 包含 以 下 儿 个 步骤 :恒星 选取 、 获 取 误 差 数 据 、 建 立 指向 误差 模型 、 


计算 模型 参数 、 


难 验证 模型 


4 符合 度 、 在 望远镜 控 和 


上 程序 中 导入 修正 系数 并 实时 修正 。 


3 常用 指向 误差 模型 


(1) 球 谐 函数 模型 。 


称 为 球 谐 函 数 模型 ， 


它 可 应 用 于 任何 机 架 结 构 的 望远镜 , 通过 
合 ， 方 法 简单 ， 可 拟 合 各 种 误差 。 对 地 平 式 望远镜 和 赤道 式 望 远 镜 ， 
度 。 但 其 模型 参数 较 多 且 没 有 实际 物理 意义 ， 


数据 后 ， 求 出 的 模型 参数 变化 大 。 
(2) 基本 参数 模型 。 对 望远镜 的 静态 指向 误差 及 基本 参数 逐 
归 因 子 。 由 于 导致 指向 误差 的 因素 较 多 且 可 能 存在 相关 性 ， 通 过 逐 


中 的 不 显著 因子 剔除 ， 
个 引入 ， 每 引入 一 个 变量 


采用 球 谐 函 数 多 项 式 对 地 基 望 远 镜 指向 误差 进行 拟 合 的 误差 模型 


步骤， 直至 添加 所 有 误差 源 。 此 时 ， 


过 逐步 回归 分 析 ， 故 参数 项 较 少 
的 机 架 结 构 引 入 其 他 参数 ， 设 计 不 同 的 模型 并 验证 残 差 ， 

(3) 机 架 模 型 。 
结构 分 析 建 立 模型 。 全 模型 共有 24 项 参 
19 项 参数 ， 其 中 方位 轴 参 数 有 11 项 ， 高 度 轴 参数 有 10 项 (方位 胃 
西向 倾斜 为 共同 参数 ) ， 这 
和 高 度 轴 各 项 参数 。 


机 架 模 型 也 和 彩 


模型 
参数 之 


确保 模型 中 的 参数 项 必须 存在 
旦 ， 要 对 模型 


J 


尔 为 转台 模型 ， 


对 基准 本 


口 


Sd 


且 不 可 替代 。 


为 球面 的 误差 进行 拟 
能 达到 较 高 的 内 符合 精 
参数 之 间 相 关 性 大 ， 模 型 不 稳定 ”"， 重 新 测 得 


步 回 归 分 析 以 确定 显著 回 
步 回 归 分 析 ， 可 将 误差 源 
其 分 析 思 路 是 将 误差 源 逐 
FPF 的 原 有 变量 逐个 检验 ,将 不 显著 的 因子 删除 ， 重复 该 
FPF 的 所 有 误差 源 对 修正 系数 都 是 显著 的 。 因 为 
辣 相 关 性 小 ， 模 型 较为 稳定 ”。 但 需 根 据 望 远 镜 


过 程 复杂 且 可 拓展 性 较 差 。 


些 参数 均 有 实际 的 物理 意义 


过 对 参数 间 的 相关 性 


， 表 1 和 表 


表 1 机 架 修 正 模型 参数 表 (方位 轴 ) 


Tab. 1 Table of parameters of frame correction model (azimuth axis) 


常用 于 地 平 式 望 远 镜 ， 通 过 对 望远镜 机 染 
器 归 分 析 及 试验 ， 保 留 
南北 向 倾斜 和 方位 轴 东 
2 所 示 分 别 为 模型 方位 轴 


序号 


误差 项 描述 


方位 函数 


1 


2 


3 


CCD 视 场 中 


方位 轴 南 北向 倾 人 


ca S 
四 以 大 


P 心 与 光 胃 


位 轴 与 高 度 轴 非 正 交 


圆 率 (sin) 


E33 


椭圆 紊 (cos) 


1 


-cosAtanH 


-sindtanH 


secH 


-tanH 


Sind 


cosd 


序号 误差 项 描述 方位 函数 
8 方位 上 的 双 周 期 (经 验 项 ) sin24 
9 方位 上 的 双 周 期 〈 经 验 项 ) cos24 
10 方位 轴 的 双 周 期 比例 (sin) sin2AsecH 
11 方位 轴 的 双 周 期 比例 (cos) cos24secH 


表 2 机 架 修 正 模型 参数 表 (高 度 轴 ) 


Tab.2 Table of parameters of frame correction model (height axis) 


序号 误差 项 描述 高 度 函 数 

1 高 度 编码 器 零点 差 1 

2 方位 轴 南 北向 倾斜 sind 

3 方位 轴 东 西向 倾斜 -cosA 

4 镜 简 弯 沉 coth 

5 高 度 编 码 器 的 阻力 (sin) sind 

6 高 度 编 码 器 的 阻力 (cos) cos4 

7 高 度 轴 的 阻力 (sin) Hsind 

8 高 度 轴 的 阻力 (cos) Acos4 

9 高 度 轴 的 双 周 期 (sin) sin24 

10 高 度 轴 的 双 周 期 (cos) cos24 
机 架 模 型 的 表达 式 为 : 
AA=ai~a2zcos Atan H ~a3sin Atan H +a4ssec H ~—astan H +assin A 

+a7COS A+assin2A+ascos2A+aiosin2Asec H+auncos2AsecH 1) 

AH =bitb2sin A—b3cos A+ba4acot H +bssin A+becosA ed 


+b77H sin A+bDbsH cos A+bosin 2A +biocos 2A 


机 架 模 型 是 对 基本 参数 模型 的 扩展 ， 具有 更 高 的 精度 ,符合 度 优 于 基本 参数 模型 。 由 于 
机 架 模型 加 入 了 更 多 的 参数 ， 增 大 了 参数 间 的 相关 性 ， 因 此 没有 基本 参数 模型 稳定 ”。 


4 误差 数据 获取 与 机 架 模 型 系数 求解 


望远镜 指向 模型 修正 分 为 零点 校正 ,指向 误差 测定 和 模型 系数 求解 三 个 模块 ， 并 将 修正 
系数 载 入 控制 程序 , 在 恒星 和 卫星 的 指向 与 跟踪 中 实时 修正 。 本 文 基于 地 基 激 光 发 射 望远镜 ， 
选用 机 架 模 型 修正 指向 误差 。 

(1) 测定 码 盘 零 点 偏差。 通过 FK5 星 表 或 星 图 软件 让 望远镜 指向 任意 一 颗 恒 星 《〈 避 免 
指向 双星 系统 , 条 件 允 许 时 选 北极 星 可 便于 跟踪 , 避免 引入 更 多 误差 ), 手动 将 该 星 调 到 CCD 


中 心 ， 用 码 盘 测 得 的 方位 /高 度 值 减 去 该 星 的 理论 方位 /高 度 值 ， 即 为 码 盘 的 零点 偏差 。 通 过 
这 一 步骤 ， 修 正 码 盘 零 点 和 坐标 系 零 点 之 间 的 偏差 。 

(2) 划分 天 区 ， 获 取 数 据 。 在 睛 夜中 ， 根 据 拟 获取 的 测试 数据 组 数 均匀 划分 天 区 ， 在 
每 个 分 区 内 选 定 一 颗 亮 度 适中 恒星 ( 太 亮 容易 曝光 过 度 , 使 获取 的 误差 数据 不 准确 ， 太 上 暗 则 
可 能 与 周围 的 恒星 混 清 )。 此 时 有 两 种 方法 获得 指向 误差 数据 ， 方 法 一 : 根据 该 星 偏 离 CCD 
中 心 的 脱 计量 和 CCD 像素 的 空间 分 辨 率 ， 计 算 望 远 镜 在 该 空间 指向 处 的 指向 误差 。 方 法 二 : 
手动 将 目标 星 调整 到 相机 成 像 中 心 , 该 方法 能 避免 方位 轴 像 素 分 辨 率 随 高 度 角 变化 带 来 的 影 
响 ， 有 效 提 高 测量 精度 。 通 过 测量 星 点 在 CCD 上 偏 移 一 定 像素 对 应 的 望远镜 转角 ， 可 以 测定 
像素 分 辨 率 。 经 试验 ， 在 高 度 角 为 0” 时， 方位 像素 分 辨 率 见 表 3: 

表 3 CCD 像素 分 辩 率 测量 表 


Tab. 3 CCD pixel resolution measurement table 


方位 高 度 
角度 1208. 0” 894. 0” 
CCD 像素 数 2592 1944 
像素 分 辩 率 0. 466” 0. 460” 
则 高 度 角 卫 处 ， 方 位 的 像素 分 辩 率 为 (0. 466XcosH“。 因 此 ， 采 用 方法 二 记录 误差 数据 


〈 方 位、 高 度 、 方 位 误差 、 高 度 误差 )。 试 验 共 测 得 66 组 数据 。 
3) 求解 指 向 模型 系数 。 利 用 Matlab 软件 对 获取 的 误差 数据 进行 处 理 , 计算 修正 系数 。 
误差 数据 包含 4 项 ， 分 别 为 方位 角 4Z、 高 度 角 已 、 方 位 偏差 9e7ta4Z1、 高 度 偏差 de7 ta 
(方位 、 高 度 以 弧度 表示 ， 偏 差 数 据 以 角度 表示 ， 共 测 得 组 数据 , 直 1,2…, 们 ， 将 其 导入 
Matlab 中 ， 构 造 四 组 一 维 列 向 量 。 代 入 机 架 模 型 的 方位 函数 与 高 度 函 数 ， 得 到 修正 系数 矩 
阵 4， 通 过 最 小 二 乘法 求解 ， 可 求 得 修正 系数 向 量 及 推导 过 程 如 下 : 
deltaAZ 
deltaH 


X 
令 系数 矩阵 4 -| , 误差 矩阵 DET -| | 为 方位 函数 矩阵 ， 了 为 高 度 函 


数 和 矩阵 ， 根 据 机 架 模 型 中 系数 的 物理 意义 ， 有 
X=[1 0 -cos(AZ).x tan(H) -sin(AZ).x tan(H) sec(H) -tan(H) sin(AZ) cos(AZ) 0 
sin(2.x AZ) cos(2.x AZ) 0 0 0 0 sin(2.x AZ).xsec(H) cos(2.x AZ).xsec(H) 0 0] 


(3) 
Y=[0 1 sin(AZ) -cos(AZ) 0 0 0 0 cot(H) 0 0 sin(AZ) cos(AZ) 
H.xsin(AZ) H.x cos(AZ) 0 0 sin(2.x AZ) cos(2.x AZ)] (4) 


采用 多 项 式 曲线 拟 合 ， 拟 合 曲线 方程 为 19 阶 多 项 式 。 测 试 数 据 共 66 组 ， 系 数 和 矩阵 4 
为 132X19， 误 差 矩 阵 DET 为 132X1。 
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19 
det1= Dax' =aAL1 + asAn2 + +aAL1s: +aAi1 
i=0 
19 | 
det > = aa = ai9A2.1 ais42 十 +a242182 + qiA2.19 
i0 (5) 


19 
det 132 = = ai9A132.1 +ais4iaz 2 十.… 十 024132 .18 + A1A132.19 
i0 


用 矩阵 形式 可 以 表示 为 : 


DET=AB (6) 
B 为 待 求解 的 修正 系数 向 量 [bio,b1i8,…,b1] ,在 方程 两 边 左 乘 A, 得到: 


A'DET=A’AB (7) 
在 新 方程 的 两 边 同 乘 44 的 道 和 矩阵 ， 有 


B=(A’A)'A’'DET (8) 
上 式 右 半 部 分 各 项 均 已 知 ， 可 求解 得 到 修正 系数 向 量 及 
在 望远镜 的 基 座 、 重 力 分 布 未 发 生 较 大 变化 的 情况 下 ， 一 次 误差 测定 可 长 期 使 用 。 


5 误差 数据 分 析 及 模型 符合 度 验证 


在 指向 模型 数据 处 理 的 过 程 中 , 修正 系数 的 拟 合 精度 受 单 次 观测 的 统计 量 影响 。 当 出 现 
目标 星 识别 错误 ， 基 座 振动 等 情况 时 ,测量 的 误差 值 会 大 幅 偏离 真实 误差 值 。 异 常 的 观测 数 
据 在 图 形 上 呈现 孤立 型 , 偏离 由 大 多 数 数据 组 成 的 分 布 趋势 , 需要 吻 除 这 些 数据 以 降低 对 拟 
合 精度 的 影响 。 
通过 残 差分 析 的 方法 剔除 偏离 过 大 的 数据 ， 从 而 保证 实验 数据 的 精确 度 。 在 上 节 的 求解 
模型 中 将 方位 误差 提取 出 来 ， 单 独 求解 ， 由 式 1 可知， 方位 修正 系数 共 11 项 。 若 要 求解 这 
11 个 系数 ， 则 最 少 需 要 11 组 方程 。 试 验 共 获 取 66 组 数据 ， 可 构建 66 组 方程 ， 相 应 的 就 有 
66 个 残 差 。 本 文中 , 第 i 组 数据 的 残 差 e 为 该 组 的 预测 值 与 观测 值 yi; 之 差 , 即 @i= 一 yi。 
代入 上 节 求 解 模型 ， 有 : 


e=|(XxB)-deltaAZ |x3600 (9) 


上 式 残 差 的 单位 为 角 秒 ， 取 残 差 95% 置 信 区 间 ， 可 得 方位 残 差 图 ， 以 此 估计 预测 的 误差 
与 实际 误差 是 否 一 致 ， 从 而 验证 所 建 模型 的 符合 度 ， 方 位 残 差 图 如 图 1 所 示 : 


残 驮 / 角 炒 


方位 残 差分 析 图 (95% 置 信 区 间 ) 
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30 40 
方位 误差 组 序 


方位 残 差分 析 图 89% 置信 区 间 ) 


残 差 / 角 和 炒 


60 


图 1 方位 残 差分 析 图 


Fig.1 azim 


| 


uth residual analysis chart 


图 注 : 左 


图 1 的 横 多 


间 的 差 值 ， 以 


交 ， 则 说 明 获 取 的 指向 误差 数据 能 够 较 好 地 符合 模型 。 反 之 ， 若 置信 


为 95% 置 信 区 间 ， 右 图 


为 99% 置 信和 区间 


30 40 
方位 误差 组 序 


E 标 表示 测试 文件 中 数据 的 序号 , 纵 坐 标 代表 残 差 , 即 拟 合 值 与 实际 测量 值 之 


秒表 示 。 可 看 出 ， 残 差 随 机 分 布 在 零点 的 两 侧 ， 若 残 差 的 置信 区 间 与 零点 相 


数据 可 视 为 异常 数据 ， 需 要 剔除 。 以 方位 残 差 95% 置 信 区 间 为 例 ，66 


常 ， 需 要 剔除 以 提高 模型 精度 ， 


为 了 验证 指向 模型 对 1 


作为 评价 指标 ， 该 指标 可 说 


误差 表达 式 为 : 


区 间 侦 离 零点 ， 则 该 组 
组 测试 数据 中 有 4 组 异 


前 差 数据 的 修正 是 否 具 有 良好 的 精确 度 , 本 文采 用 均 方 误差 (1152) 
F 价 数据 的 变化 程度 ，N5E 越 小 ， 说 明 模 型 


述 数 据 越 准确 。 均 方 


N 


1 2 
MSE = 一 1 一 六 
Ne 


为 了 在 数量 级 上 更 直观 ， 可 对 均 方 误差 


i=1 


RMSE = NMSE 


下方， 即 均 方 根 误差 (RMSE)， 有 : 


(10) 


(11) 


以 方位 的 指向 误差 为 例 ， 分 别 计算 在 不 同 组 数 下 的 均 方 根 误差 ， 如 表 4 所 示 : 
表 4 不 同 数量 测试 组 下 计算 的 均 方 根 误差 (方位 角 ) 


Tab. 4 Root mean square error calculated under different number of test 


groups (azimuth) 


测试 数据 组 数 均 方 根 误差 ( 角 秒 ) 测试 数据 组 数 均 方 根 误差 ( 角 秒 ) 
12 737966.61 25 25.56 
13 1044325.26 30 29.41 
14 98457.58 35 30.06 
15 1253.46 40 32.45 
16 406.68 45 34.65 


测试 数据 组 数 均 方 根 误差 ( 角 秒 ) 测试 数据 组 数 均 方 根 误差 ( 角 秒 ) 


17 171.53 50 43.86 
18 19.08 55 46.60 
19 19.88 60 47.01 
20 21.78 66 46.76 


由 于 拟 合 方位 误差 需要 确定 11 个 系数 ， 故 最 少 需要 11 组 测试 数据 。 从 表 4 可 看 出 , 在 
一 定 范围 内 ， 增 加 测试 组 别 能 有 效 提高 模型 的 符合 度 。 当 测试 组 数 超过 20 组 时 ， 拟 合 数 据 
的 准确 度 反 而 随 测试 组 数 的 增加 而 降低 。 因 此 , 不 经 处 理 的 增加 测试 组 数 并 不 能 使 拟 合 精度 
得 到 优化 。 添 加 新 的 误差 数据 后 须 经 残 差 分 析 ， 吻 除 偏离 较 大 的 值 ， 才 能 降低 拟 合 误差 。 

在 已 获取 指向 误差 数据 的 情况 下 ， 可 通过 残 差分 析 ， 将 残 差 95% 置 信 区 间 与 零点 不 相交 
的 数据 依次 剔除 ， 如 图 2 所 示 。 并 计算 剔除 后 数据 的 均 方 根 误差 ， 如 表 5、 表 6 所 示 。 
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图 2 剔除 部 分 偏离 较 大 数据 后 的 残 差 分 析 图 〈 方 位 


Nt 


Fig.2 residual analysis chart after removing some value with large deviation (azimuth) 


HU 
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图 注 : 图 a 为 66 组 数据 图 b 为 62 组 数据 
取 95% 置 信 区 间 。 


c 为 58 组 数据 ， 图 d 为 54 组 数据 。 均 


表 5 有 剔除 残 差 偏离 数据 后 的 均 方 根 误 差 〈 方 位 角 ) 


Tab.5 Root mean square error after eliminating residual deviation data (azimuth) 


测试 数据 组 数 均 方 根 误差 ( 角 秒 ) 测试 数据 组 数 均 方 根 误差 ( 角 秒 ) 


66 46.76 58 26.25 


62 34.35 54 22.42 


由 式 2 可 知 ， 高 度 修正 系数 共有 10 项 ， 故 拟 合 高 度 误差 至 少 需要 10 组 数据 。 同 理 ， 对 
高 度 方 向 做 残 差分 析 。 


数据 组 
中 


0 5 
残 差 / 角 秒 


图 3 高 度 角 残 差 数据 分 布 


Fig.3 height residual data distribution 


1 


组 数 ， 由 


注 : 图 为 残 差分 布 直 方 图 ， 横 坐标 表示 残 差 ， 纵 坐标 表示 某 一 残 差 范围 的 数 拓 


< 


IHh 


图 可 看 出 ， 高 度 残 差 近似 服从 正 态 分 布 。 


高 度 残 差分 析 图 (95% 置 售 区 间 ) 


高 度 残 差分 析 图 (35% 置 信和 区间) 
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图 4 剔除 部 分 偏离 较 大 数据 后 的 残 差分 析 图 《高 度 角 ) 


Fig.4 residual analysis chart after removing some value with large deviation (height) 
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图 注 : 图 a 为 66 组 数据 ， 图 b 为 62 组 数据 ， 图 c 为 58 组 数据 ， 图 d 为 54 组 数据 。 均 


取 95% 置 信 区 间 。 


表 6 有 剔除 残 差 偏离 数据 后 的 均 方 根 误差 〈 高 度 角 ) 


Tab. 6 Root mean square error after eliminating residual deviation data (height) 


测试 数据 组 数 均 方 根 误差 ( 角 秒 ) 测试 数据 组 数 均 方 根 误差 ( 角 秒 ) 
66 5. 41 58 3;37 
62 4. 29 54 2:77 


使 用 乌鲁木齐 的 地 基 激 光 发 射 望远镜 试验 , 在 0” 一 360” 范 围 内 选择 60 颗 恒 星 ， 测量 
指向 误差 数据 ， 删 除 残 差 侦 离 较 大 的 数据 点 后 ， 计 算 的 两 轴 均 方 根 误差 如 表 8、 表 9 所 示 ， 
相应 残 差 分 析 图 如 图 6、 图 7 所 示 : 
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图 5 恒星 指向 数据 的 A-0-E 坐标 分 布 图 


Fig.5 A-0-E coordinate distribution of star pointing data 
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X 间 ) 方位 残 差分 析 图 (95% 置 信 区 间 ) 
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图 6 剔除 部 分 偏离 较 大 数据 后 的 残 差 分 析 图 〈 方 位 角 ) 


Fig.6 residual analysis chart after removing some value with large deviation (azimuth) 


图 注 : 左 图 为 60 组 数据 ， 右 图 为 56 组 数据 。 


高 度 残 差分 析 图 85% 置 信 区 间 ) 高 度 残 差分 析 图 95% 置 信 区 间 ) 
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图 7 剔除 部 分 偏离 较 大 数据 后 的 残 差分 析 图 《高 度 角 ) 


Fig.7 residual analysis chart after removing some value with large deviation (height) 


图 注 : 左 图 为 60 组 数据 ， 右 图 为 56 组 数据 。 


表 7 应 用 机 架 模型 计算 的 剔除 残 差 偏离 较 大 数据 前 后 的 修正 系数 对 比 表 


Tab.7 Comparison table of correction coefficient calculated by frame model (before 


and after removing data with large residual deviation) 


序号 数据 优化 前 数据 优化 后 
方位 高 度 方位 高 度 

1 0. 19809 0. 04877 0. 19549 0. 05176 

和 0. 03339 -0. 02918 0. 03242 -0. 02667 

3 -0. 00849 0. 00125 -0. 00649 0. 00380 

4 -0. 13096 -0. 01822 -0. 13219 -0. 02299 

5 -0. 00400 -0. 02918 -0. 00641 -0. 02667 

6 -0. 00100 -0. 00125 0. 00395 -0. 00380 

7 -0. 00076 0. 01993 -0. 00280 0. 01674 

8 0. 00352 -0. 00593 0. 00628 -0. 00099 

9 -0. 00143 -0. 00241 0. 00022 -0. 00224 
10 -0. 00218 -0. 00457 -0. 00349 -0. 00332 
11 0. 00007 -0. 00029 

表 8 剔除 残 差 偏离 数据 后 的 均 方 根 误差 方位 角 ) 

Tab.8 Root mean square error after eliminating residual deviation data(azimuth) 

测试 数据 组 数 均 方 根 误差 ( 角 秒 ) 测试 数据 组 数 均 方 根 误差 ( 角 秒 ) 
60 9.46 56 7.79 
表 9 剔除 残 差 偏离 数据 后 的 均 方 根 误 差 〈 高 度 角 ) 
Tab.9 Root mean square error after eliminating residual deviation data(height) 

测试 数据 组 数 均 方 根 误差 ( 角 秒 ) 测试 数据 组 数 均 方 根 误差 ( 角 秒 ) 


60 18.47 56 11.39 
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试验 获取 了 lm 人 卫 激 光 测 距 望 远 镜 , 地 基 激 光 发 射 望远镜 等 项 目的 多 组 指向 误差 数据 ， 
通过 计算 分 析 得 出 结论 : 剔除 残 差 偏 离 的 测试 组 能 够 有 效 提高 指向 修正 数据 的 拟 合 精度 。 

在 外 界 环 境 未 发 生 明 显 变化 的 前 提 下 , 指向 误差 数据 可 分 次 获取 (可 将 分 批 获取 的 数据 
合并 ， 使 用 残 差 图 分 析 ， 若 残 差 集 中 在 95% 置 信 区 间 内 且 未 出 现 断 层 ， 则 可 判定 外 部 环境 未 
发 生 明 显 变化 )， 添 加 新 的 测试 数据 时 ， 要 保证 添加 后 测试 组 的 均 方 根 误差 小 于 原 有 测试 组 
的 均 方 根 误 着 ， 方 能 提高 指向 模型 对 误差 数据 描述 的 准确 度 。 
通过 多 元 线性 回归 检验 机 架 指向 修 正 模 型 对 偏差 数据 拟 合 的 精确 度 ， 如 表 10 所 示 : 

表 10 机 架 模 型 检验 表 


Tab. 10 Frame model inspection form 


让 


测试 组 数 R”(R 为 相关 系数 ) F 统 计量 F 检 验 的 p 值 误差 方差 
66( 方 位 角 ) 0.99695889 1803.05339 2.0198 X10% 2.0253 X10™ 
54 (方位 角 ) 0.99891837 3971.18443 2.3422x10-60 4.8698 X10” 
66〔 高 度 角 ) 0.99951715 17151.7968 1.0660x10-93 2.5792X10” 
54 高度 角 ) 0.99986543 48827.2601 8.4612x10”7 6.9740X107 


R 为 决定 系数 ， 其 值 越 大 ， 拟 合 程度 越 高 。 由 表 10 可 看 出 ， 原 始 的 66 组 测试 数据 拟 合 
优 度 超过 99%， 误 差 拟 合 精度 非常 好 。 通 过 残 差 分 析 筛 选 数据 后 ， 拟 合 优 度 又 得 到 了 进一步 
提高 。 同 样 的 ， 也 可 通过 检验 的 p 值 验证 指向 模型 的 有 效 性 。 一 般 的 ， 当 p<0. 001 时 ， 说 
明 模型 有 效 ， 由 上 表 可 看 出 ,对 于 所 列 的 测试 组 , 检验 的 p 值 远 远 小 于 0. 001, 说 明 指向 模 
型 能 够 很 好 的 拟 合 偏差 数据 。 
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